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ABSTRACT 


Various computational methods and operational computer codes used to predict the 
aerodynamic coefficients and separation trajectones of aircraft stores are examined. The 
semi-empirical aeroprediction code Missle DATCOM is used to obtain the coefficients 
of a modeled store. These coefficients, together with the modeled ejection forces, are 
used in free-stream state-space equations of motion to predict the store trajectory. The 
results are compared with the Nielson Engineering and Research (NEAR) store separation 
code which provides accurate trajectory profiles, for speeds below the critical speed, by 
use Of a vortex-lattice and panel method. Modification of the Missile DATCOM 
aerodynamic coefficients provides single-point state-space prediction of the store pitch 
trajectory within 30% of the NEAR code results. Store trajectories were restricted to the 


first 0.2 seconds of free flight. 
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i STORE SEPARATION INTRODUCTION 


A. BACKGROUND 

The prediction of the trajectory of a store ejected from 
an aircraft is of major concern to defense aviation. The 
increasing requirement for aircraft to fulfill a multiple role 
in the hostile environment demands that it also carry an ever- 
increasing variety of stores. These stores range from 
missiles and bombs’ stockpiled for many years to new, 
aerodynamically complex stores with lifting bodies. 

The aim of this work is to review current methods of 
predicting store releases, compare results of two different 
methods, and to present an easy-to-use methodology for 
investigating the subsonic release of a unsophisticated store. 

A non-axisymmetric pod ejected from the F/A-18 outboard 
pylon is presented as an example. This store was chosen 
because it 1s of practical value due to the extensive use of 
these modified pods at the Pacific Missile Test Center, at 
Point Mugu, California. Non-axisymmetrical stores are very 
difficult to model and present a special problem in the store 
separation field. This pod is, therefore, modeled as 
axisymmetric and non-axisymmetric, and the results are 
compared. In addition, the store lacks aerodynamic control 


surfaces and is highly unstable. Application of the 


methodology presented here should eventually provide the user 


with accurate results for a minimal time and cost expenditure. 


B. MATHEMATICAL MODELING 

There are many ways to mathematically model. store 
separation. Economics and accuracy are always the main 
considerations when exploring the codes and methods to be 
used. For simple cases, intermediate approaches such as panel 
methods and solutions of the Euler equations’ provide 
sufficient accuracy. Calculations involving flow separation 
and shock wave interference become too complex (for 
intermediate methods and acceptable results can only be 
obtained using modern computational fluid dynamic (CFD) 
techniques involving solutions of the Navier-Stokes equations 
(Ref. l:p. 7-2). For speeds below the subsonic Mach critical 
speed, component buildup methods (empirical and semi- 
empirical), and panel methods provide the necessary accuracy 
and are much easier to use. The codes chosen for this study 
are explained in detail in their respective sections. 

The F/A-18 model was obtained from Mr. L. L. Gleason of 
the Ordnance Systems Department, Naval Weapons Center, at 
China Lake, California. The F/A-18 was chosen due to its wide 
variety of store carriage and also because of its solid future 
with Naval aviation. 

The store model is based upon both the Missile Datcom 


method [Ref 2] and also by the Nielson program [Ref 3]. These 


will be discussed in Sections III and IV. Appendix A contains 
the physical description of the pod and its inertial 


characteristics. 


C. METHODOLOGY 

The initial investigation (Section II) consisted of 
looking at the ejector forces, combined with the store’s 
inertial characteristics, to give a sense of the magnitude of 
initial velocities and moments. These forces and moments, 
while easy to predict, were insufficient to provide real 
evidence of a safe release or of a release problem. A large 
pitch-down moment or a large vertical velocity may be off-set 
by unforeseen aerodynamic forces. The velocities and moments 
calculated will be used in Section III and compared with 
results of Section IV for trajectory prediction. 

Section III introduces a method of predicting the store’s 
longitudinal and lateral aerodynamic coefficients using a 
missile code (Missile DATCOM) developed by McDonnell Douglas 
Missile Systems Company for the Flight Dynamics Laboratory 
(FDL), Wright-Patterson Air Force Base (WFAFB), Ohio [Ref. 2]. 
These aerodynamic coefficients are then used in longitudinal 
and lateral equations of motion to predict the free-stream 
trajectory of the store. A brief explanation of the 
derivation of these equations of motion is provided to the 
reader. An accurate representation of the aerodynamic 


characteristics is essential for the full prediction of the 


shape’s trajectory. Missile DATCOM is based on the body 


buildup method and includes a number of prediction methods for 


each component of the configuration. Other codes, such as 
MISL3, are also available, each with its own relative 
strengths. 


Section IV describes a computer prediction method 
developed by Nielson Engineering and Research, Inc., (NEAR), 
also under contract to the United States Air Force. This code 
provides a six degree-of-freedom (6DOF) simulation which takes 
into account the aircraft flowfield using vortex-lattice and 
panel methods. The NEAR simulation provides a high fidelity 
representation for the subsonic case, but also requires the 
most effort and computer capability. An attempt has been made 
to minimize the complexity of the input to the program. Once 
a good data base has been developed the NEAR code should prove 
easy to use. It has been used extensively throughout the 
defense industry, and has undergone many modifications to 
incorporate improvements and options. 

Section V compares the trajectory results of the 
linearized aerodynamic simulation method and the NEAR code. 
A modification to the linear coefficients used in the linear 


aerodynamic simulation is also discussed. 


II. EJECTION FORCES AND MOMENTS 


Although quite uncomplicated in nature, calculating the 
ejection forces and moments provides the user with a feel for 
the magnitude of the initial movement of the store. This 
procedure should be done prior to using any other method as a 


preliminary step in order to prepare for the simulation. 


A. EJECTOR CHARACTERISTICS 

The non-aerodynamic forces involved in the ejection of a 
store depend upon the ejector cartridge, the ejector rack, the 
weight of the store, and the rigidity of the wing. Each of 
these parameters determine the resultant moments and forces 
provided to the store. The store and the wing (or fuselage) 
were considered as rigid, thus neglecting any aeroelastic or 
structural bending effects. None of the methods discussed 
here address aeroelastic bending due to the added complexity 
of tne problem. 

The Douglas production BRU-32/A Bomb Ejector Rack combines 
two sets of hooks, one set at 14-inch spacing and one set at 
30-inch spacing, with an ejection system designed for carriage 
of stores with suspension lugs per MIL-A-8591. Two 
electrically initiated CCU-45B cartridges are used for store 


release and ejection. The self-retracting ejector pistons, 


Spaced symmetrically 20 inches apart at each end of the rack, 
have a piston stroke of six inches. These pistons are spring 
loaded against the store during loading to prevent impact of 
pistons during firing. The orifice sizes can be varied, by 
replacement, to provide force and pitch control for store 
separation. [Ref. 4] 

The best source of ejector data is from the Aircraft 
Ordnance Procedures (AOP) contained in the Aircraft Stores 
Interface Manual, (Reference 4). This manual contains a 
complete description of the rack under consideration. Here, 
the example separation uses two 14-inch spaced hooks and 
Q0.118-Diameter orifices. The force diagram corresponding to 
this orifice, modeled from the AOP, is shown in Figure 1. 
Sufficient data points were read from the AOP graph in order 
tee curve-fit the points using Computer Associates’ 
CricketGraph® software [{Ref. 5]. This provided two things; 
1.) integration of the polynomial equation gave the total 
impulse value, and 2.) Nielson software in Section IV uses a 
fifth-order polynomial for calculations. Note that the AOP 
shows the total force provided by the ejector, while the NEAR 
Simulation requires a force-per-ejector-foot polynomial. 
Given the shape of the force curve, a triangular 
representation of the curve gave a value of 348 lb,-sec versus 
the integrated value of 344.9 lb,-sec, which corresponds to an 
error of only 1%. Therefore elementary estimates of the total 


impulse are practical. 


Ee es _ 
7000 a /o 

6000 

5000 

4000 


3000 3 - 


FORCE (LBF) 


2000 eo i 
1000 


0 -« pone = 
0.000 0.010 9.020 0.030 0.040 0.050 
TIME (SECS) 


Figure 1 Ejector Force 


B. FORCES 
1. Theory 

Using the conservation of linear momentum theory 
presented in Appendix B, the pod was modeled as a simple beam. 
The weight was concentrated at the center of gravity and the 
resultant ejection forces were placed at the ejector feet 
locations. The ejection force time histories are treated as 
an equivalent force impulse which results in a change in 
linear momentum as shown by the following equation: 

JIF, dt = G, - G, = (M*V,), -(M*V,), 

The integral for representing the linear impulse is 
based upon the fifth-order polynomial force time history shown 
here. 
y= 280.2-608022.8x + 180688545.8x* - 10638451765.5x° + 

234526098234.3x* - 1783873760363. 4x? 
2. Example 

For our example, the BRU-32A bomb rack used two CCU- 
45B cartridge-activated devices (CADS). The peak force was 
14,500 lb, and the pulse duration was 50 msec. The total 
linear impulse calculated was 345 1lb,-sec. Using a weight of 
371 pounds for the pod, the end-of-stroke velocity was 29.9 
feet per sec (fps). This value was used in Section III as the 


initial velocity for the trajectory simulation. 


The graphs contained in the AOP shows end-of-stroke 
velocity for a 371 1b store to be approximately 24.5 fps. 
This graph is derived from empirical data and should provide 
the most accurate representation of the actual release 
velocity. The values of 24.5 and 29.9 fps are compared with 
the end-of-stroke store velocity of 30.9 fps predicted from 


the NEAR code in Section IV. 


C. MOMENTS 


1. Theory 
The pod was modeled as a simple beam. The inertial 
characteristics of the store are listed in Appendix A. The 
two ejector feet provide the beam with a moment. The 


equations used are elementary, and with certain assumptions 
the angular velocity can be calculated. The equations are 
based on conservation of angular momentum. The equations and 
assumptions are presented in Appendix B for completeness. The 


resulting equation is shown here. 


Pee eer eal 

2. Example 
The center of gravity was offset aft of the center of 
the ejector feet by 4.05 inches. The impulse thus provided 
the pod with 1.12 radians per second (rps) or 64.1 deg/sec 
initial angular velocity. The question now is, are these 
reasonable values? at the end of 0.429 seconds, the shape 


will be one body-length below the aircraft. At the same time, 


the shape will have rotated 27.5 degrees nosedown. Therefore, 
without considering the aerodynamics of the vehicle, the pod 
seems to have cleared the wing. A good rule of thumb is 2-3 


body lengths clearance, although exceptions do occur. 


D. PRELIMINARY OBSERVATION 

This initial calculation provides an approximate estimate 
of the forces and moments involved. The values obtained are 
of sufficient accuracy for calculations such as end-of-stroke 
loading, etc. Excessively large separation velocities are 
usually an indication of a miscalculation. Translational 
velocities should range from 10-30 fps and pitching velocities 


from U0z0 "hp s- 
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III. AERODYNAMIC CALCULATION 


A. THEORY 

Many physical systems can be modeled by second-order 
differential equations. The mathematical treatment of fixed- 
wing flight vehicle motions was first developed by G.H. Bryan. 
He laid the mathematical foundation for airplane dynamic 
stability analysis, developed the concept of the aerodynamic 
stability derivative, and recognized that the equations of 
motion could be separated into a symmetric longitudinal motion 
and an unsymmetric lateral motion. Experimental studies were 
initiated by lL. Bairstow and B.M. Jones of the National 
Physical Laboratory in England, and by Jerome Hunsaker of the 
Massachusetts Institute of Technology to determine estimates 
of the aerodynamic stability derivatives in Bryan’s theory. 
In addition to determining stability derivatives from wind- 
tunnel tests of scale models, Bairstow and Jones 
nondimensionalized the equations of motion and showed that, 
with certain assumptions, there were two independent 
solutions, i.e., one longitudinal and one lateral. These two 
solutions provide the free-stream trajectories we seek. [Ref. 
op. 113] 

The dynamic stability characteristics of a pod can be 


represented by six equations of motion, three for the forces 


At 


involved X, Y, and Z@ and three for the moments L, M and N. 
The force equations relate the forces acting on the body to 
the corresponding linear accelerations and the moment 
equations relate the moments to the corresponding angular 
accelerations. It is usually possible to consider the 
longitudinal motions completely separately from the lateral- 
directional motion, by neglecting the various coupling terms. 
(Ref. 7p. 14] 

As a caveat to the use of this method, these equations are 
the linearized version and are only valid up to approximately 
TOredegrees anglesotederdek. Treatment of the non-linear 
aerodynamic coefficients, while not extremely difficult, does 
require knowledge of the behavior of the coefficients. Since 
the goal of this report is to predict the behavior of new 
shape configurations, this knowledge is not presumed. 
Therefore, the results obtained in this fashion are only valid 
for estimating the motion in the first fractions of a second. 
These results are compared with those obtained in the more 
rigorous method of Section IV. 

1. Longitudinal Equations 

The rigid-body, longitudinal equations of motion can 
be developed from Newton’s second law: 
L2Pitching moment s==M,,=1,,*© 
The assumption that body motion consists of small deviations 


from its equilibrium flight condition allows us to use 
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perturbation theory to examine the aerodynamic force 
derivatives in terms of angle-of-attack (AOA), vertical 
velocity, pitch angle, and pitch rate, by means of a Taylor 
series expansion. The X-force, Z-force, and pitching moment 
equations comprise the longitudinal equations. To separate 
these equations from the lateral equations, they must not be 
coupled. This is a reasonable assumption given the nominal 
geometric shape of missiles or pods [Ref. 6]. The body is 
constrained to move in a vertical plane and is free to pitch 
about its center of gravity. For a comprehensive derivation 
of the equations of motion, see Reference 6. The resulting 
equations of motion are shown in Appendix C. 
2. Lateral Equations 

The Y-force, rolling, and yawing moment equations 
comprise the lateral equations of motion. Once again, these 
equations are uncoupled from the longitudinal equations by the 


assumption of small cross-component moments of inertia. That 


erel.,, L:;, ty. acre small in comparison to the principal axis 
moments. These equations are also valid assuming only small 
variations in displacement and velocity. The lateral 


equations are derived from the following Newton’s laws: 
2Rolling moment s=1,,*® 
LYawing moments=I,,*y 

The third equation is derived from the Taylor series expansion 


of the side force derivative. The lateral motion of the pod 


es 


disturbed from its equilibrium state is a complicated 
combination of rolling, yawing, and sideslipping motion 
Assuming the cross products of inertia are ignored, some of 
the coupling terms can be simplified. 
3. State Variable Representation 
The linearized longitudinal and lateral equations 
developed above are simple, ordinary linear differential 
equations with constant coefficients. The coefficients in the 
differential equations are made up of the aerodynamic 
stability derivatives, mass, and inertia characteristics of 
the pod. These equations can be written as a set of first- 
order differential equations in the state-space (state 


variable) form: 


{x} = [A]*{x} + [B]*{n} 


where {x} 1s the state vector, {N} 1s the control vector and 
the matrices [A] and [B] contain the pod’s dimensional 
stability derivatives. In our case the pod does not have any 
active control surfaces. For missile launches, however, the 
control vector and the B-matrix would be used to represent the 
control surfaces and control input. The simplified state- 
space representations of the longitudinal and lateral 


equations of motion are shown in Appendix C. The required 
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coefficients listed in Appendix C are derived using the 
procedures in Section III.B. 

The state-space representation of the equations of 
motion can be solved simultaneously using matrix software such 
as MATLAB® [Ref. 8] or Control-C® [Ref. 9]. Given the 
conditions of the ejection, namely, initial translational 
velocities, initial angular velocities, and attitude angles, 
solutions of the state-space equations can be used to predict 
the trajectory of the pod. The initial angles and velocities 
are the initial conditions imposed upon the state-space 
equations. The Control-C® commands for this procedure are 
contained in Appendix D. 

B. MISDATCOM 
MISDATCOM was developed by McDonnell Douglas Astronautics 
Company, St Louis, Missouri, for the Flight Dynamics 
Laboratory of the Air Force Wright Aeronautical Laboratories, 
Wright-Patterson AFB. The program was completed in December 
poe 5 . The Missile Datcom was created to provide an 
aerodynamic design tool which has the predictive accuracy 
Suitable for preliminary design, and the capability for the 
user to eaSily substitute methods to fit specific 
applications. [Ref. 2:p. 1] 
1. Theory 
There are many different methods to predict the 


aerodynamic static, dynamic, and control characteristics of 


10s: 


missiles. Component build-up was chosen as the most suitable 
for this program. Although panel methods are better suited 
for arbitrary configurations, component build-up was chosen 
due to the accuracy provided for conventional configurations 
and for the ability to do parametric studies easier. 
Basically, component build-up consists of using 
various methods to compute the characteristics (skin friction, 
force and moment coefficients, panel loading,...) of the 
individual configuration components. The various methods are 
chosen for their applicability to the configuration or flight 
Gendieren:. Then the components are combined. Previous 
methods of combining the components (fins, body, engine inlet) 
consisted of adding the individual coefficients and then 
multiplying the sum by some interference factor obtained using 
slender body theory. The approach taken with MISDATCOM was to 
use the "equivalent angle of attack method" developed by 
Nielson Engineering and Research, Inc. (NEAR). This method 
assumes that the panel loading for a given panel angle-of- 
attack is unique. With this method the panel angle of attack 
is computed including the effect of panel roll orientation 
with respect to the free stream velocity vector, panel 
proximity to the fuselage or to other panels, and external 
vortex flow field effects. Then the isolated panel 
characteristics are interpolated at the panel equivalent angle 
of attack to yield the panel load when mounted on a body in 


combination with other surfaces. [Ref. 2] 
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2. Procedure 

The procedure for using MISDATCOM is straightforward. 
The pod shape is modeled using simple geometric shapes. The 
previous release of MISDATCOM software could only handle 
axisymmetric shapes or forms with a vertical plane of 
symmetry. Unfortunately, the pod under consideration is non- 
axisymmetric and therefore the results are not entirely valid. 
The latest version, however, can model some non-axisymmetry 
through the use of a "protuberance" option. This is the April 
1991 release and is now available. 

For this investigation, a Simulation run was made 
using both versions and a comparison of the aerodynamic 
coefficients was made. The pod radius was also varied from 
the minimum pod to the maximum, and the resulting coefficients 
compared. The MISDATCOM code coefficients were somewhat 
insensitive to radius changes of this magnitude, therefore, 
the minimum-radius, axisymmetric case was retained for 
comparison with the minimum-radius non-axisymmetric case. 
Figure 2 shows the pod semi-profile of the MISDATCOM procedure 
alongside the pod semi-profile of the NEAR simulation. Figure 
3 contains a sample input to the program and Figure 4 lists a 
partial sample output. The actual input and output of the 
program are contained in Appendix E. These aerodynamic 
coefficients are next used in the longitudinal and lateral 
state-space equations of motion to predict the store 


trajectory. 
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Figure 2 MISDATCOM Pod Shape 
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RENs4.I13E6, 
MALPHA-&. ,ALPHALO. Of) ogee 4, BIIhe-0.5,8 
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PRINT AEKO HOUY 
PRINT GEOM BONY 
NEXT Case 


A WARNING A DHE KERERKENCE AREA IL UNSPECIFIED, HEEAULT VAIUL ALEUMEN 
& UVARNING & THE REFERENCE LENGTH IS UNSPECIETED, OFEFAULT VALNE ASSUMED 
THE BOUNDAFY LAYER IC ASCUMFD 10 BF TUKRULENT OVER ALL CUMPONEHITS OF THE CONE IGUPATION 


THE INPUT WNITS ARE IN INCHES, THE SCALE FACTOR Is 1.0000 


Figure 3 MISDATCOM Input 
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C. AERODYNAMIC CALCULATION 
1. Longitudinal 

The aerodynamic coefficients derived from the 
MISDATCOM were entered into the plant matrix [A] of the 
longitudinal equations of motion. The initial conditions were 
applied and time history of the pitch angle, pitch angular 
velocity, and vertical velocity was found. Figures 5-8 show 
that the pod assumes a large nose-down, highly divergent 
pitching motion. This is understandable since there are no 
control surfaces to make the pod stable. Examination of the 
roots of the plant matrix [A] indicates an unstable flight 
vehicle. The pitch angle and pitch rate shown are in the 
store body coordinate frame. 

These results are valid only for the range of linear 
values of the aerodynamic coefficients. This is approximately 
up to 10-12 degrees angle-of-attack. Therefore, an estimate 
of the non-linear behavior of the pod after an AOA of 10 
degrees is reached is necessary. Due to the large pitch 
angular velocity, it is obvious that this limit is reached 
after only 0.2 seconds. Section V discusses the acceptable 
time range where these results are valid. 

To estimate the non-linear behavior, the dominating 
terms in the [{A]-matrix must be determined. This provides an 
insight into the possible range of values to substitute into 


the [A]-matrix. A discussion of a possible approach to this 
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Figure 5 Pitch Angle Response 
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Figure 7 Vertical Separation Distance 


VERTICAL VELOCITY (FPS) 


~®-Pod 2 27d 
-@-Pod 1 2d 


4 #05 O08 O07 O08 09 10 
TIME (SECS) 


0.4 





Figure 8 Vertical Velocity 
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problem is presented in Section V. To completely investigate 
this area is beyond the scope of the method outlined here. 
Comparison of the "linear" trajectory obtained here will be 
made with the trajectory obtained in Section IV, which does 
take into account the non-linearity of the coefficients. 
2. Lateral 

Due to the nearly vertical forces and moments provided 
by the ejection rack during straight and level flight, the 
free-stream lateral equations of motion will not provide us 
with any insight into the safe jettison of the. shape. 
However, 1f the pod were experiencing sideslip, then the 
sideslip could be entered into the state-space equations as an 
initial condition. For this pod a sideslip of -0.5 degrees 
was assumed. The resulting lateral motion is not’ shown 
because the main emphasis is on the longitudinal motion. The 
pod is unstable laterally, also, but the initial movement is 
small due to the relatively small initial conditions. In 
addition to the yaw angle and rate, if the ejection rack 
provided an initial roll rate, such as when the pod is loaded 
off-center of the bomb rack, then that influence could also be 


included. 


D. PARAMETER VARIATION 
The parameter variation due to non-axisymmetry is 


difficult to handle. Comparison with shapes with known 


Ze 


(experimentally obtained) coefficients might provide some 


degree of accuracy in predicting more accurate results. 


E. PRELIMINARY OBSERVATION 

The results of the method in Section II were a vertical 
velocity of 29.9 feet/sec and an angular velocity of 64.1 
deg/sec. The aerodynamic method of Section III begins with a 
vertical velocity of 29.9 fps, quickly diverging to an 
extreme value. The pitch rate, gq, also diverges qQuickERe 
Because the pod is very unstable any linear approximation of 
its behavior will have strict limits. The coefficients used 
in this calculation, from the MISDATCOM code, are not valid 
beyond a fraction of the trajectory. The valid time range of 


prediction is presented in Section V. 
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re "NEAR" PROGRAM 


A. BACKGROUND 


A computer prediction method was developed by Nielsen 


Engineering and Research, Inc., (NEAR) under contract to the 
United States Air Force. The work was performed during the 
Perrod 1968 to 1972. The final result is a method for 


predicting the six degree-of-freedom store separation 
trajectory at speeds up to the subsonic critical Mach number. 
After delivery of the program to the Government, many new 
capabilities were added. The code used for this paper was 
obtained from the Naval Weapons Center, China Lake, Ca. The 
program has been widely accepted by industry and government. 
The code encompasses 9,900 lines of code and thus requires a 
device with sufficient computer memory for operation. (This 
also depends upon the program application.) 

The aircraft fuselage, separated store body, and adjacent 
stores are modeled using point sources and sinks. Angle of 
attack effects are included using a cross-flow model. The 
aircraft wing and wing pylons are modeled using planar vortex 
lattice models which include dihedral, camber, and twist of 
the aircraft wing. Thickness strips are used to model the 


thickness of the aircraft wing and pylons. [Ref. 10:p. 807] 


Z5 


The capability exists to install multiple sets of wings, 
fins, or canards, and to use active control surfaces on the 
store by inputting the control laws into the program. Powered 
separations may be simulated by inputting the store thrust 
Characterigercs 

The NEAR program actually consists of two separate 
programs: the source program and the trajectory program. Both 
are described in Section IV.B below. 

Alternate separation programs include USTORE and USSAERO 
codes. USSAERO was developed by F. A. Woodward, of NASA, as 
a lower-order panel method. USTORE was developed from USSAERO 
by G. J. van den Broek, of the National Institute for 
Aeronautics and Systems Technology, Pretoria, South Africa. 


[Refs Lip. S09) 


Be THEORY 

The three principal tasks in the prediction of a store 
trajectory are: first, the determination of the nonuniform 
flow field in the neighborhood of the ejected store; second, 
the determination of the forces and moments on the store in 
this flow field; and third, the integration of the equations 
of motion to determine the store trajectory [Ref. 3]. 

1. Source Program 

The source program 1s used to represent an 

axisymmetric body as a distribution of sources along the axis 


of the body. It provides point source-sink distributions to 
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represent the fuselage, rack, and store volumes. The program 
calculates and prints the source strengths and locations. 
These quantities are then used as input data to the trajectory 
program. 

The source program is used for the generation of an 
aircraft or pylon model which is then used for the trajectory 
program. These models are Mach number dependent and thus need 
to be generated for each test case with different Mach number. 
For ongoing store separation studies, a good database of 
aircraft models is required and should be available from 
appropriate Government facilities. The F/A-18 model was 
obtained from China Lake along with the program code. This 
model included pylon stations. 

2. Trajectory Program 

The trajectory program uses the source-sink 
distributions from the source program, and additional 
information to first determine the vorticity distribution 
which represents the wing-pylon loading including interference 
of tne fuselage, rack, and stores, and then to calculate the 
store trajectory. See Sections IV.C.3.d.(1) and (2) for 
descriptions of the vortex lattice method and the panel 
method. 

Once the trajectory program input has been generated 
for a single flight condition, it is relatively easy for the 


user to input different stores. The major effort is to obtain 
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the initial aircraft fuselage, wing, and pylon models for this 


input. 


3. Shape Modeling 
a. Aircraft 

The aircraft fuselage geometry was modeled by using 
a Fortran program called NGDELX developed by L. Gleason of the 
Naval Weapons Center. This program provides the coefficients 
to the NEAR polynomial representation of equivalent body. The 
shape is divided into appropriate segments. Radius values for 
the segment points and maximum radius values are entered. 
These coefficients are then entered into the source program to 
generate the source representation of the fuselage. Figure 9 
shows the coefficients for the F/A-18 fuselage. Following is 


the equation used to calculate the coefficients: 


Coan = Cy + C,* (C,* (K/L)* 4 Cy*(X/L) +e, + Cy* (X/L) + Co) 


b. Pylons 


The pylons used are modeled the same as in Section 


Iv.C.l1.a above. The geometry is used to obtain the NEAR 


coefficients. 
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Polynomial Representation of Equivalent Body 


-.01456 





Figure 9 F/A18 NEAR Polynomials 


c. Pod or Missiis 
Again, the procedure in Section IV.C.1.ais used to 
represent the shape as a set of polynomial coefficients. The 
NGDELX program is again used to represent the missile shape as 
an "equivalent body of revolution" (EBR). NGDELX is very easy 
to use, but, again, only represents axisymmetric shapes. 
d. Wing (Lifting $Svusteces) 
The aircraft wings, missile fins, or any other 


lifting surfaces are modeled using a vortex lattice method to 
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represent the loaded wing. The thickness is also modeled 
using the lattice method. 

(1) Vortex Lattice Methods. There are several 
variations of the vortex lattice method that are presently 
available and have proven to be very practical and versatile 
theoretical tools for the aerodynamic analysis and design of 
planar and non-planar configurations. [Ref. 12 p. 27] The 
vortex lattice method represents the wing as a planar surface 
on which a grid of horseshoe vortices is superimposed. The 
velocities induced by each horseshoe vortex at a specified 
control point are calculated using the law of Biot-Savart. A 
summation is performed for all control points on the wing to 
produce a set of linear algebraic equations for the horseshoe 
vortex strengths that satisfy the boundary condition of el 
flow through the wing. The vortex strengths are related to 
the circulation and the pressure differential between the 
upper and lower wing surfaces. The pressure differentials are 
integrated to yield the total forces and moments. [Ref. 13 
De 61} For a rigorous introduction to the vortex lattice 
method, see Reference 13. 

(2) Panel Methods. The configuration is 
modeled by a large number of elementary quadrilateral panels 
lying either on the actual aircraft surface, or on some mean 


surface, or on a combination thereof. To each elementary 


panel, there is attached one or more types of singularity 
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distributions, such aS sources , vortices, and doublets. 
These singularities are determined by specifying some 
functional variation across the panel (e.g., constant, linear, 
quadratic, etc.), whose actual value is set by corresponding 
strength parameters. These strength parameters are determined 
by solving the appropriate boundary condition equations. Once 
the singularity strengths have been determined, the velocity 
field and the pressure field can be computed. [Ref. 13:p. 258- 


Za 9 } 


C. PROCEDURE 
1. Source Program 
a. Input 

The program input consists of the polynomial 
representation of the equivalent body of revolution (EBR) 
values obtained in Section IV.B.3.a. The surface of the EBR 
1s approximated by these polynomials, which represent the EBR 
x,r distribution. The source program is then run with a user- 
specified finite distribution. The surface obtained from this 
source distribution can then be compared with the polynomial 
surface. Additional runs may be required to closely match the 
surfaces. The aircraft fuselage, pylon(s), and the adjacent 
store EBR values are input along with two variables, NRAT and 
PERCR. NRAT is the number of segments which the body will be 


divided for the specification of the source distribution. 


Sil. 


PERCR is the source spacing for each NRAT segment and is input 
as a fraction of the local body radius of the segment. 
Adjacent stores are included in the aircraft source 
LnPUEs In this way it becomes part of the aircraft 
configuration. The separated store is not input to the source 
program since the trajectory program calls the source program 
during its execution to model the store. this is due to the 
fact the store changes position during the program run and the 
sources/sinks will change values. Figure 10 lists a sample 


format for the source program input. 


FILE: 2-Fissou T 


1f-18 AIRCRAFT 
l 
1 
F-18 AIRCRAFT 
6 0.70 
o00000000000 00000000 
-12 0.28 48 -542 
35763 -~1.0 -292589 ~12790 
-01978 0.0 90 -0 
-03239 0.0 .0 .0 
-88523 -1.0 -82195 
-99071 =1.0 28 
0 ~0.01456 .018683 
001 1.1 ~ 6400 
6 
- 10046 0.15525 0.28311 0.37443 0.97717 
0.6 0.8 1.0 1.2 8 


0 
0 
0 
0 
0 
l 
0 
0 
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Figure 10 Source Program Input 


S2 


5b. Output 
The source program output is used as input to the 
trajectory program. However, it is not directly read into the 
trajectory input and must be entered via keyboard. The form 
of the output is shown as a partial output in Figure 11. 


Source locations and strengths are listed. 
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P-10 AIRCRAFT 


X/L OF END POINT OF BACH SECTION OF 8007 


SECTION 1 2 $ 4 § ¢ 


xu! 0.12000 0.20000 0.48000 6.34700 


6.44000 1.00000 


COEFFICIENTS OF POLYNOMIAL DESCRISINO EACH SECTION 


SECTION cl c2 cs 4 cs cé Cc? 
i -0.35768 -1.00000 0.29259 0.12790 0.00000 0.00000 1.00000 
2 0.01970 6.00000 6.00000 0.00000 0.06745 0.00000 0.900000 
8 0.03239 0.00000 0.00000 0.00000 60.023te 0.00000 6.00000 
4q 0.08528 -1.00000 0.02198 0.34430 0.00000 0.00000 -1.00000 
0 -0.99071 -1.00000 1.29000 0.40496 0.00000 @.00000 1.00000 
4 0.00000 -0.01¢65¢6 0.01063 -0.00280 0.00000 0.00000° 1.00000 


FIRST SOURCE AT X/L= 0.00100 

LAST SOURCE AT M/L © 1.10000 

FROM 6.00100 Tc 6.10066 @OURCE SPACINO 18 0.60000 TIHES LOCAL RADIUS 
PROM @.100¢¢4 To 9.19328 S0URCR sPactnHa 19 @.90000 TIHES LOCAL RmaDtuUsS 
FROM @.15328 To 0.20311 SOURCE SPACING 13 1.00000 TIHES LOCAL MaDIUS 
FROM @.20311 TO @.2974¢43 SOURCE SPACIHO 18 1.20000 TIHES LOCAL AADIUS 
PROM @.37443 Ta O.97717 G@OURCE SPACING 18 0.90000 TIHES LOCAL AADIUS 
FROM O.97717 TO 1.00000 ¢0UACE SPACING 15 ©.60000 TIHES LOCAL AADIUS 
FOR THIS CASE THERE ARE 74 SOURCES 


Figure 11 


SOURCE LOCATIONS AND @ODY RADIUS AND 


aL 
A/L 
DR/DX 


x/L 
AsL 
DR/DX 


M/L 
R/L 
DR/ DX 


R/L 
AsL 
CR/DX 


K/L 
AsL 
DR/OK 


RAL 
R/L 
OR/OX 


R/L 
RZL 
DAs OX 


RAL 
R/U 
DR/OM 


RAL 


0.00100 
0.00041 
0.40501 


0.00308 
6.00125 
6.39906 


@.0093¢ 
6.00368 
0.37904 


0.02716 
6.00994 
O.32418 


@.07072 
6.02130 
0.19946 


@.15993 
@.03087 
0.06746 


O.36176 
0.0403) 
0.023180 


0.82264 
0.06008 
0.13330 


O.73649 


0.00118 
0.000<8 
0.40524 


0.003861 
@.00146 
@.39734 


@.01092 
0.004726 
0.37406 


@.03142 
0.01129 
O.SL189 


6.07980 
6.02301 
@.1748? 


0.10176 
6.03204 
@.06748 


0.37630 
@.04i11 
6.02318 


@.$S011 
@.0S168 
0.00624 


@.76749 


0.00138 
0.00056 
6.40487 


0.00424 
6.00171 
0.39533 


0.01275 
@.0069¢ 
@.3¢6030 


0.03625 
@.01277 
@.29709 


6.00971 
0.02460 
@.14605 


0.20464 
0.033390 
@.06745 


@.39979 
6.04166 
@.02310 


@.87961 
0.08376 
0.08702 


0.79794 


0.00162 
0.00066 
0.40379 


0.00497 
0.00199 
0.39290 


0.01406 
0.003973 
O.3¢6167 


0.04172 
0.01435 
@.20012 


0.10025 
0.02601 
0.12003 


6.22063 
0.033920 
0.06745 


0.42359 
@.04221 
6.02310 


0.61032 
0.05800 
0.02030 


0.02771 
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34 


0.00190 
0.00077 
0.40267 


6.00588 
@.00233 
0.39024 


6.01733 
6.00661 
6.35408 


0.04707 
@.01¢02 
0.2634) 


@.11139 
@.0z2719 
0.09069 


0.25377 
6.03690 
0.06745 


0.44771 
0.04277 
6.023180 


O.64n7s 
0.05549 


-0.00088 


0.05665 


SURFACE SLOPE AT THESE LOCATIONS 


0.00224 
0.00091 
0.40100 


0.00602 
@.00272 
0.30705 


0.02017 
@.007¢0 
0.345394 


0.05474 
O.01776 
0.24390 


6.12693 
@.029 34 
@.06745 


@.20012 
6.03090 
@.02310 


0.47214 
0.04333 
0.023180 


0.67349 
0.05534 


-0.00007 


0.00660 


0.00268 
0.00106 
6.40084 


0.00799 
6.003t7 
0.38334 


0.02842 
0.000871 
0.335406 


6.06235 
0.01994 
6.22257 


0.14312 
0.02948 
@.06745 


0.30708 
0.03953 
6.02310 


0.49690 
@.04508 
0.1022? 


0.70511 
0.05493 
-0.01732 


0.91144 


2. Trajectory Program 

The trajectory program uses the output of the source 
program, and the input of the wing configuration and pylons, 
along with store information to calculate the trajectory of 
the separated store. It is not the intent of this report to 
fully explain the aircraft modeling details of the program. 
The assumption is made that the appropriate aircraft model is 
amamlable for thescorrect flight condition and configuration. 
The applicable store inputs to the program are discussed and 
the resulting trajectory examined. Nielson contains a 
thorough description of the theory and input of the aircraft 

model in Reference 3. 

a. Input 

The example input to the trajectory program was 
1064 lines in size. The store data entry begins at 
approximately 980 lines, depending upon aircraft input. The 
store input will vary according to the physical configuration 
of the store. For the example separation, Input Items 35-43 
and 48-51 were used to input store physical data. The ejector 
information was input in Items 44A-47. This is where the 
fifth-order polynomial representation of the ejector force was 
used. Data concerning initial and final times, and 
integration step sizes were entered in Input Item 72. The 
correct input form of the above data is contained in the NEAR 


code users manual [Ref. 14] and is not presented here. 
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b. Output 

The output of the trajectory program is dependent 
upon the input conditions, 1.e., the number of stores, number 
of store segments, and number of integration steps, etc. The 
example case output is 4,524 lines. The first 1,840 lines of 
output are aircraft-specific output, such as the source and 
VOrEITCICY IMNfLOrmaelon, 

The last 2,650 lines of code contain the trajectory 
data. The ejector force data 18 repeated as is the store 
shape data, reference dimensions and inertial characteristics. 

Proceeding the ejector data, a data block is 
presented for each time integration step. This data block 
contains the parameters which describe the store displacement, 
velocity, and acceleration during the separation. A partial 
sample output of this section of the output file is included 
mm  Fagure = ize Plotting routines can be developed to 
facilitate data analysis of the simulation but have not been 
developed for this-~ study. The following are short 
descriptions of some of these parameters. A comparison with 
the results of the method of Section III follows in Section 
ny. 

The force and moment coefficients are listed by 
their individual contributions. The effects of buoyancy, 
slender body theory, crossflow, and empennage are totaled to 


give an effective C,, Cy, Cy,, and C,,. 
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8000 BJECIOR DATA 


eroot NEroLY(t? KECET Towra (ed grace ery 
a 4 -t.0008 @.90 9e?. 0008 
Trewoqe tt TPENOC ETD TPPHIOCE. Sf eRe. 4p TPEIOVI 9D 
‘ ®.eegeee-618 @.,e0e0e ee @,¢0¢08-°e@@ © ,e¢¢e0@¢>o@ @.¢066F°-6e@ 
eoccee TEHE 10 THE INDEPENDENT BF VECTOR VARIABLE evcece 
Foor wNerourvsel) weete tee taeity Steowe rt?) 
e 1 @.ese @.0@ 899.9008 
Teepe. gy) TRESIHCLE.Ot Teerme«s.9t TPeRNOIt. a) TEPVIDGE. St 
@.eeece-@) @.ecoce ee @.e0cce-ee @.CCCoe-ce @.e0cee -ce 
eecee TIME JO THE INDEPENDENT BJECTOR VARIABLE eeuce 
COSCFFICIENTO OF OTH DEGREE POL YHOMINAL 
€FooT a ct ce ce cq ce ce 
8 1 @.20098-09 -@, ¢0eor-es @.tesvecer ~@,.foeeeeyn @.29688°02 -@.)7e4R°48 
2 2 O.2R 028-09 ~0. den 08°06 @.10078°0F -@.J00¢e+g) @.9tase-s2 ~0.37048°49 
CTORE NUMBER 1 89 THe Grone feces 
ADDITECHIAL BHFUT POR INIG Qrore 
@TORNE HAZE @ BE.$06 ¢LU09 
MOMENTS AND PRODUTIS Of BHEAITEA, ILVO - 80 PT 
ix e 9.00 
byYy e $2.00 
{37 @ $2.06 me 
YZ e @.08 
XZ © @.0@ 
ixy e@ 6.00 : 
STORE MOMENT CENIER 19 -8.400 FEET RENEID NOSE 
OLORE AEFEREIICE AREA 19 2.499 Se.Ft, 
O@TOME AEFEACHCE LEWOTH 89 1.499 ef. 
STONE CENTER OF DAAVITY OFPSET EROM MOMENT CENTER, FEET 
ROAR @ §.000¢8 
YOAR e 8.000086 
7BAR © @, 00008 
POL YIOHIALS SFECIF YEP COMPRESOTELE @fone gare 
M/L OF EHD OF FACH GECTION 
eection wl TtvFe PRY 10 e@ HO. B ce POLYT 
9 0.01990 e 
2 @.08900 8 
9 0.097¢0 Q 
4 0.13¢00@ 0 
3 @.79470 0 
e 1.900098 e 
COEFF ECRENTO CP POLYNOMIALS DESCAL@ING EACH SECTION 
eeCcTion el c2 cs ce cs ce c? 
1 @.ecece 6.90008 @.eecce @.cccce @.¢decce @.ccece @.09000 
2 @.0t219 -8.009000 @.rreee -@.e0rTIs @.cccee @.eccee 1.80008 
8 -@.00179 -§.00000 0.99900 -@.00¢0f2 9.00000 ©. 00000 1.90000 
4 0.01048 -4.00008 O.27010 -0.01940 0.00000 @,90000 1.00000 


Figure 12 Trajectory Program Fartial Output 
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TIME « @.tf30 Ther e 1.0000 line e @.0908 
Teg © ©, 19300 sEcorios 


PORCE AHD MORTEHNT COEFFICTONT® TREF, AREA © 2.090 Pee? -— AEF, LEMOIN «© §.633 FI.9 
cn cy c.M cLNn cul 


OvUOVAIRY ~0 s0007 ©.01039 ©.09019 -0.081926 
SLENDER @00Y 1.76770 -@.02¢079 @.9eler -@. 10729 
CROSSFLOW @.10708 -@.00243 -0 49002 @.0cee?7 
CHF Ereiwo€ 3 6.00000 0.00008 0.60008 @.00008 @.600006 
TOTAL 1.97906 -@.00061 @.39900 -0.22969 6.00000 


LOCATION OF @TORE BI FUSELAGE COOADIHATE @YOTEN. DIMENSIONS OF FEET 


RELATIVE TO FPUSELACE NOSE RELATIVE TO UHTTTAL FOSTT ION 

uP vr zP CFL KF DEL VF DEL 2° 

HOS® ~27. 96876 -10. 097 OF @.00908 -0.00228 -@. 00107 4.19358 
MoMOM ~889.86922 ~-11.09016 @.19066 -@.00029 -@.000186 8.32506 
@OASE -60.60001 -11. 10601 0.997237 @.0eses @.00T19 O.G7071 


TRANSLATIONAL VELOCITIES AND ACCELERATIONS OF GIORE IN FUSELACE COORDINATE OveTEM 
RELALIVE TO FUSELADE MOTION 
rxP CYP oP DIP O2vP O22° 
0.16904 ~-8.00000 UF.07662 ~-1.20077 -6.00393 21.099278 


ROTATIONAL VELOCETIZE® AND ACCELERALIONG OF @TORE 1M STORE COORDINATE OYSTIEN 7 
P e a Poot pool root : 
@.000278 -$.01096 -0.00260 6.00000 0.64030 -@.00762 


OTORE ANOULAR OR TONTATION IN PUSELACE CDORDIHATO OYOTEN AHO AATO® OF CHAHOC® OF THESE AOL Oe 
AJIOLOG 1M DOGDCEES. AATOO OF CHIANO® [NH RADIANS FER OECD 
Pet THETA Pe or st OMA orn 
-08.01029 -@.71236 0.002766 -@.00279 -3.01096 @.00069 


TivMe «@ @.0390 TIHEP © 1.0000 €TiKe « @.0900 
TIME © 6.19008 8ecoNDeE 


PORCE AND HOMEMT COEFFICIENT®O IREF, AREA © @.099 PYeeg -- REP. LENOIN eo 3.638 P7.3 


CN cv CLR CUN cuUL 
QUOYAIKY ~@.0ce7@ @.00708 @.10979 -08.G4708 
SLEHOER GODY 0.74699 -0.02490 @.090270 -@.07090 . 
CROSSFLOW 0.16919 ~@.00220 -8.60s00 @.ccete 
CHP EHNADE | 6.00008 6.00000 @.00008 @.00000 @.00000 
TOTAL 1.36309 -@. 00763 @.27310 -0.210386 @,00008 


LOCATION OF BIORE IN PUSELADE COORDINATE OYOLOM, DIMENSIONS OF FEET ¥ 


ROLATIVE 10 PUSELADE HOSE ROLATIVE 10 INTTIAL POSIT Tor 

uP vP 2 ot Cel wR bel @F 

WOSE ~-27.8069F -05.08063 B.e9sse ~ -0.10990 ~-@.00141 4.9470 
HOM 6-8 3.26077 -11.19000 0.04094 «9 @.0T400 -8.00072 @.u9sra 
OASE ~40.07729 -h0. 10069 4.47720 0.10666 @.00150 8.19668 


TRANSLATIONAL VELOCETIE® AND ACCELERATIONG OF OTORE IM FPUSELACE COORDINAIE OvSIEH 
AELATIVE TO PUSELADE MOTION 
OxP Over OIF O2xP O2vP 02 2P 
0.00099 ~-@©.00199 U9.01869 ~§. 270760 -0,00309 S8f.05107 


ROTATIONAL VELOCEITIC® AND ACCELERATIONS OF BIORE IN OTORE COORDINATE SYOTEM 
° 8 R Foo! poor AooT 
@.00090 -@.98703 ~0.00739 @.00000 0.69669 -@.00693 


OTORE ANOULAR OR TENTATION IN PUSELAGE COORDINATO SYSTEM AND RATER OF CHAHOE OF TISE ANCLOO 
ANOLOG@ 16 DEDROLO. RATO® OF CHANDOS IH RADIANS PER SECOCHHD 
Pes THOTA PHI or et OMetA dP NI 
-O.00360 -00.0929) O.0030% -@.00318 -@.38700 0.00003 


TOME ec @.9950 7IMEP «© 1.0008 @7IH@ «© 6.0908 
Tita © @.07900 eecoros 


PORCE AND MOMENT COEFPICIENT® UREP, ARZA © @.090 FIsef ~~ REP, LENOIN © 0.988 PT.1 
cM cv CLN CUM Cu 
GUO YANCY ~0.000¢0 @.00799 @.10089 -0,00468 
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The ejector force and moment components are listed. 
The vertical force history matches the force profile from 
Section II very closely. The store’s nose, inertial reference 
center, and the base position in the fuselage coordinate 
system is listed. The separation distance from the initial 
position is also listed as output. Plots of these parameters 
will indicate the miss distance of the shape from the wing. 
The store translational and rotational velocities are output 
in both the store’s reference frame and the fuselage reference 


frame. 
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Ne, RESULTS AND COMPARISONS 


A. EJECTOR FORCES AND MOMENTS 
1. Vertical Velocity 
The predicted end-of-stroke vertical velocity from the 
hand-calculation was 29.9 fps, comparing with the NEAR code 
results of 30.84. This is a 3% difference. The predicted 
value from the ejector graph of Reference 4 yields 24.5 fps. 
The ejector graph is representative of the average store used 
in the Naval inventory, which more than likely would have 
control surfaces for stability reasons. These surfaces would 
also provide aerodynamic damping which could be the source of 
the discrepancy. 
2. Pitch Rate Q 
The predicted pitch rate, q, of 1.12 radians/second 
matches well with the NEAR’ end-of-stroke q of 1.08 
radians/second. For most applications, this correlation of 4% 


should have sufficient accuracy. 


B. TRAJECTORY COMPARISON 

The main consideration in determining the accuracy of 
these trajectory methods is the safe separation (or jettison) 
of the pod. Vertical distance and pitch attitude are 
indications of the pod separation while velocity and angular 


velocity are indications of the store loads. Therefore, 
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Criteria for judging the correlation between the linear method 
and the NEAR code must be based on the practical standpoint of 
a reasonable assumption of safe separation. The individual 
parameters are addressed below. For determining the valid 
time range of the aerodynamic prediction method, the NEAR code 
was used as reference. 
1. Vertical Separation Z 

Vertical distance z in the linear method is derived 
from the vertical velocity term. Therefore the error 
(deviation from NEAR results) will be slightly less due to the 
integration effects. However, given the initial vertical 
velocity, the NEAR code predicts a five foot separation. The 
error limit must be set as an acceptable variation of this 
distance. Based upon practical knowledge of store separation, 
a 20% error limit was set. Figure 14 compares the linear 
aerodynamic method results comparison with the NEAR results in 
the amplified time region of the first 0.2 seconds. 

2. Pitch Angle 

Since the pitch angle is used solely for determining 
the safe separation and is not used for more accurate 
applications, like rocket engine ignition, the effect of an 
error in pitch will vary as the sine of the angle deviation. 
This could result in an error, due to rotation about the 
center of gravity, of an order of 7.2 feet times the sine of 


the error angle. The pitch angle error must be fairly small 
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initially, but can grow as time increases due to increased 
separation distance and its decreasing effect. Therefore, a 
conservative error limit is set for 0.1 sec. At 0.1 seconds, 
the NEAR pod has rotated 4.5 degrees. A 30% error would give 
a vertical error of 0.2 feet, which 1s approximately 10% of 
the actual separation. Figure 15 shows a comparison of the 
linear method with the NEAR results. Figure 16 provides the 


error between the methods as a function of time. 
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Figure 14 Vertical Separation Distance 
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Figure 15 Pitch Angle A 
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3. Pitch Rate 
The, “secuinred .AaGcuUracy for “piteh rate, will be 
determined by its application. Pormeunictance, navigational 
considerations will require quite accurate predictions of the 
angular loading of its components. However, for this case, 
there is not a requirement of this kind. The divergent nature 
of the pitch rate from the linear aerodynamic model does not 
present a problem for this case, beyond the fact that the 
pitch angle follows this trend. Therefore, an error limit is 
not placed upon the pitch rate. 
4. Vertical Velocity 
Vertical velocity is divergent and represents the lack 
of aerodynamic damping in the linear model. Therefore, it 
does not provide meaningful information beyond the phase 
relationship with the vertical separation distance. For this 


reason, an error limit is not placed upon the velocity. 


C. AUGMENTED AEROPREDICTION 

Due to the limited valid time limits on the aerodynamic 
method, an alternative approach was investigated. Obviously, 
the non-linear operating region is reached very quickly. The 
two dominant terms of the A-matrix are the damping terms, Cy, 
ana Cy. The unstable pod is shown from the NEAR code to 
quickly reach a somewhat steady-state value of pitch rate as 
it tumbles. An effort was made to calculate a new (and valid) 


value of G,. This approach modeled the tumbling store as a 
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cylinder of constant radius. Using a value of C,=1.0, and the 
end-of-stroke pitch rate of 1.12 radians/second, new values of 
damping coefficients were calculated. 
1. C,,- Coefficient 
The new pitch damping coefficient of -.1054 was used 
versus the MISDATCOM value of -.04. This is an obvious 
improvement and also realistic since the pod will definitely 
experience aerodynamic damping, due to pitch rate, that would 
not have been predicted by MISDATCOM. 
2. Cy- Coefficient 
The value of the aerodynamic damping effects of change 
in angle of attack, G,, was combined with C,, in the MISDATCOM 
output. Therefore, the valid value of CG, is unclear. For 
this model of a rotating pod, it was assumed that Cy, would 
have minimal effect on the damping and a value of 0.0 was 
assigned. 
3. Trajectory Prediction 
The predicted parameters using the above coefficients 
are contained in Figures 17, 18 and 19. As shown in the 
figures, there is considerable improvement in the correlation 
with the NEAR trajectories for pitch angle and pitch rate. 
The more realistic value for GQ, obviously improves the 
results. The time range of validity for these two parameters 
1s now estimated to be 1.0 seconds. However, the most 


important parameter, vertical separation distance, is still 
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outside of reasonable accuracy limits. Therefore, the 
augmented plant matrix method 1s still unusable for complete 
trajectory prediction. Figure 20 shows the parameter error 
for the augmented system. 

Although there is substantial correspondence in the 
partial comparison between the NEAR results and these results, 
there are not enough data points to substantiate this method. 


Further work using different shapes need to be done. 
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Figure 17 Vertical Separation Distance (Augmented) 
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Figure 18 Pitch Angle ‘Agomented’ 
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Figure 19 Pitch Rate (Augmented) 
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Figure 20 Parameter Error (Augmented) 
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VI. CONCLUSIONS 


In the preliminary design application, MISDATCOM provides 
an easy-to-use method to predict moderately accurate store 
aerodynamic characteristics. These aerodynamic coefficients 
can be used to model the store dynamics for design work. 
However, their use in the linearized equations of motion for 
an externally ejected store trajectory prediction is highly 
limited. 

Comparisons with the Nielson Engineering and Research 
(NEAR) code show that trajectory values are outside of 
reasonable accuracy limits for use in store separations. The 
store separation vertical distance and pitch angle values 
diverged within Q.1 seconds. Therefore, the use of the 
linearized equations of motion cannot provide useful 
trajectory information. 

However, a modification to these coefficients to more 
accurately reflect the existing aerodynamic damping does 
provide correlations with the NEAR code of within 30% for 
pitch angle and 10% for pitch rate. Use of these "augmented" 
linear equations of motion for attitude prediction shows 
sufficient accuracy for many applications. 

The NEAR code is widely used for ejected store trajectory 


prediction in the aviation field. Its use for safe separation 


0 


investigation should be considered as standard for obtaining 
accurate results. The cost of the Nielson code in both 
engineering time and computational time increases with the 
level of information produced, but the need for. safe 
separation accuracy, aerodynamic load distributions, and 
vortex-induced effects make the extra effort justified. 

For limited applications requiring less accuracy, the 
augmented linear equations of motion can provide some 
information for the initial trajectory movement. A potential 
application would be the investigation of the effects of small 


weight or inertial modifications to existing stores. 
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APPENDIX A 


A. POD CHARACTERISTICS 

The pod in question was an ALE-29B pod modified to 
accommodate the SUU-53 chaff dispenser and associated 
equipment. The ALE-29 pod was flight approved by the Aircraft 
Configuration Control Board (ACCB) document #75-004. It was 
originally an ALE-2 pod. No significant aerodynamic changes 
have been made to the pod. 


The pod physical characteristics are listed in Table 1. 


TABLE 1 POD CHARACTERISTICS 
Length 153.8 inches 


Radius —-minimum 9.75 inches 


—-maximum 11.20 inches 


CG location -empty 62.4 inches from nose 
=u LL 65.4 inches 
S20) lbs 
Sil bs 
of Inertias 

-Ixx SCOmIb— ft 

-Iyy 1675 lb-ft? 

-Izz 1675 lb-ft?’ 
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APPENDIX B 


A. FORGE 2OUATIONS 


ZF = Mtv = d/dt (Mtv) = G 


Separating variables: 


DE aC LF, = G, Dee) 
Assuming ZF, , = 0, 
Integrating the remaining equation, 
ice CMG oc MEV). = (MeV). 


Assuming straight and level flight, V,, = 0 yields: 


JIF, dt = M*Av, = MF¥V, 


Integrating the ejector force polynomial yields a total 


impulse of 345 lb,-secs. Substituting (wt=371 lb), 


345 lb;y-cecs = 371 1b, * V, 


We ipeyy a e(lpe—cee/ ib) * 32.174 (1b.-ft/1b,-sec’*) 


V 


zi Zora ps 


SS 


Be MOMENT EQUATIONS 
Mo = exe en M*v = d/dt(r x Mtv) = H. 


where H, 1s the angular momentum about pt O. Integrating the 


moments, 
JIM, dt = H+ rx H 


Separating the component equations yields, 


2M, = TRO 
XM, = I,*@, - (I, - I,,) *W,*0, 
EM ee Ore et Ler eee es 


Assuming the reference frame coincides with the principal 


axeS, Ixy xz,yz2 = 0, amd assuming 0, , = 0, yields the folUc wii 


Jim, dt = I * @, 


@, = 1.12 rad/sec 


54 


APPENDIX C 


A. LONGITUDINAL EQUATIONS OF MOTION 

The equations of motion for the pod can be derived from 
Newton’s Second Law of motion, which states that the summation 
of all external forces acting on a body must be equal to the 
time rate of change of the momentum of the body, and the 
summation of the external moments acting on a body must be 
equal to the time rate of change of the moment of momentum 
(angular momentum). The time rates of change are all taken 


with respect to body coordinates space. [(Ref.8] 


center 


Homent 





Figure 21 Store Reference Frame 


Do 


> AF,=m= (U+W*Q-V*R) 
) AF =m* (W+V* P-U*Q) 
» AM=0*1 7 Pt Gee Pe eeu 


XZ 


> AF, =m* (U+W*Q-V*R) 
\* AF,=m* (W+V*P-UxQ) 
> AM=0+I, 


By restricting the disturbances to small perturbations 
about the equilibrium condition, the product of the variations 
will be small in comparison with the variations and can be 
neglected, and the small angle assumptions can be made 
relative to the angles between the equilibrium and disturbed 
axes. 

Yo AF =m*u 


>) AF, =m* (w-U*q) =m* (w-U+®) 
\> AM=I,*q=I, 48 


Sxpanding the applied forces and moments in terms of the 
changes in the aerodynamic and gravitation forces and moments 


using the total differential form yields: 
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> dF,= (8F,/dU) *dU+ (OF,/dW) «dw+ (GF,/dW) *dW+ (OF,/00) *dB- (dF, /, 


> dF,= (OF ,/dU) *dU+ (GF ,/dW) *dw- (dF,/dW) *dw+ (AF,/08) *d@- (dF,/08 


SD A M= (GM/QU) *dU+ (OM/dW) *dW+ (OM/dW) «dwWw+ (dM/d®@) «de 


Non-dimensionalizing the terms in the above equations 


yield the 


aerodynamic coefficients. Following are the 


definitions of the longitudinal stability coefficients and 


derivatives: 


ort 


Cc 


xa 


Variation of drag and thrust with u. 

Lift and drag variations along the x-axis. 
Gravity 

Downwash lag on drag. 

BEtece Of PlECh “rate on draq. 

Variation of normal force with u. 

Slope of the normal force curve. 

Downwash lag on lift of tail. 

Effect of pitch rate on lift. 

Beeects Of thrust, slipstream, and flexibility. 
Static longitudinal stability. 

Downwash lag on moment. 


Damping in pitch. 


Arranging the linearized longitudinal equations into 


matrix form yields the following matrix equation. The control 


Sa 


matrix B iS not shown here but is included in the Control-C 
program. Since there are no active control surfaces, the B- 


Matrix 1S a ZeEro Matrix. 


Au X, xX, 0 =o Au 
Aw Z ibe u 0 Aw 
BG | | ME EM Ze reel ae Ur ee 0 Aq 
LAO | ih 0 0 0 0 A@ 
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Be LATERAL EQUATIONS OF MOTION 
In the Same manner used to obtain the longitudinal 
equations, the lateral equations are derived. The same 
assumptions using perturbation theory apply. Following are 
the equations in their developing order. 
\> AF, =m* ( V+U*R-W* P) 


) AL=P*T,.-R*J,, +O*R* (1,-1,) -P*O*J,, 
>) AN=R*I,+P*J,,+P*O* (I,-I,) +O*R+J,, 


\) AFy=m* (v+U*r+uxr) 
)) AL=p+1,-r*J,, 
Y) AN=L*I,-p*J,, 


> AF =m* (v+U*r) 
\> AL=p*I,-£*J,, 
) AN=r*I,-ptd,, 


> dF,= (OF ,/0B) *dB+ (OF, /OV) *d¥+ (dF, /dD) +d@+ (dF, /d@) «db-+ (dF, /: 


AB | Y3/u ao ~— (1-Y,/u) g*cos@8/u AB 
Ap Lg ice L. 0 Ap 
Ar| = Np N, N. 0 AY 
Ad | 0 1 0 0 | [Ad] 


De 


APPENDIX D 


m=input(’entec m (1lbs):')/32.2) 

Voeinput(’enter Uo (fps):'); 

Seinput(’enter S (in*42)2:'°)7/144; 

beinput(’entec b (in):’)/12; 

Rowinput(’enter Ro (slugs/ft**3):'); 

crefeinput(’enter c ref(in):’)/12; 

Ixeinput(’enter Ixx (ft-lb):')/32.2; 

Ty-input(’enter IYY (ft-lb):')/32.2; 

Izeinput(’enter Izz (ft-1b):')/32.2; 

g-32.174, 

Cn=-input(’enter Cn (Normal Force coefficient):’); 

Cm=input(’enter Cm (Pitching Moment Coefficient):’); 

Caeinput(’enter Ca (Axial force Coefficient):’); 

Cyeinput(’enter Cy (Side Force Coefficient):’); 

CLNeinput(’enter CLN (Yawing Moment Coeff):’); 

CLLeinput(’enter CLL (Rolling Moment Coeff):’),; 

CNAeinput(’entec CNA (Normal Force Deriv wrt AOA):');3 

CMA-input(’entec CMA (Pitching Moment Deriv wrt AOA):‘); 

Si ile | ues CYB (Side Force Deriv wrt Beta):’)} 

CLNBeinput(’enter CLNB (yawing Moment deriv wrt Beta);’), 

CLLBe einput(’ entec CLUB (Rolling Moment deriv wrt Beta):’); 

CLeinput(’enter Cl (Lift coef£):’),; 

CD-input(’enter Cd (Drag coef£):'),; 

CNQeinput(’enter CNQ (Nocmal Force Derciv wrt Q):')} 

ChQwinput(’ enter CMQ (Pitching Moment Deriv wert Q):’)) 

CNAOADOTe=input(’enter CN AOAdot (Normal Force Deriv wrt AOAdot):’); 

CMAOADOT@input(’enter CM AOAdot CES Mom Deriv wet AOAdot):'); 

q@=.5*ro*uo*uo) 

Xuw-2*Cd*q*6/(m*uo) 

Keacl Gua ttmetoh 

Zue-24*Cl *q*sa/(m*uo) 4 

Zwe-Cd*q*s/(m*uo) } 

Zwd=Cnaoadot*cref*s*g/(2*u04u04m) 3 

MweCma*q*s*cref/(Iy*uo); 

Mu=0; 

Ma-uo*Mw} 

Mwd=Cmaoadot*cref*q*atcref/(24uo*uo0*ly)) 

Mad-uo *Mwd } 

Mq=Cmq*cref*q*s*cref/(2*tuo*ly) 

Ae@=(Xu Xw 0 -g3Zu Zw uo 073 (Muthw 
(MqtMwd4u0)/(1-Zwd) 0;0 0 1 0 

B-(0 0 0 0;0 0 0 Oj}’scH{0 10 0 

1ce{0 0 0 0)’, 

IC(2)#input(’enter initial vertical velocity (fps):')) 

1C(3)einput( ‘enter initial pitching rate (rad/s):')3 

ttt-input(’enter end time (6ec)’), 

t-0: .00L:ttt, 

uuue(ttt/.oolpel, 

u-0*onese(2,uuu)) 

simu(’1IC’,IC), 

yesimu(a,b,c,d,u,t)) 

erase} 

plot(t,y): 


Zu) /(1-Zwd) (Mw+Mwd*Zw)/(1-Zwd)... 


d* 
)3 
3 


D-(0 0); 





Figure 22 Control-C Program For EOM Simulation 


60 


APPENDIX E 


0000°t S) yOLDVi JIVDS JH “SIHONI MI 3B¥v SLING INdNI FHI 


NOITIYVUNOIANOD AKL JO SINANOAWND TI¥ AZANO ANAINAANI 326 UL MA¥WNSsVY SI AZAVI sTVANNOE TKI 
GZUNSSY 3BNIWA Tinvs3ad *AALITIIASNN SI HIONTT JZINIAIIZA AHI ¥ SuINsvn ¥ 
UIWNSSY ANIwH LINwsI4d “AAILATIZASWN SI VANY JINIATAAA AWL ¥ ONINEVA ¥ 


3SWVD IX3N 
l1aqadg@ wOxD Inl ud 
; ANON OUAXIVY INIEd 
é ating 
; s°O°telivad’s°sveliyvit 
**OselIWI SSL sTINADE SO COUAINAIT* SGT eASONMS “97 SASONT AONT XS 
$°S2°S9eydx O13948 
BC * Oe VLE Oh 6°37 0* CORUNA TY YUN TUN 
; “S93CT°reNnay 
*Z2°OrPHDUN’ CLesiduWN NODITIs 
duvg 
NI wId 
ova AISIG 


IV1 TINI NOTLVANDTINAD AMON TIIWIS UIIsvo 


3S¥D SIH 803 INdNI SOBVD ZHI 3sY ONIAOTIOI 


SINJNT 3S) 
SNOLIVAMN) ING JVISSIN Sud SOOKIIN JIWYNLOORIY 


¥ 69/L AIA ¥ NUIJLYON WIISSIM AILwWOLNY AYSN JHI 


* 


MISDATCOM INPUT/OUTPUT 


Figure 23 


61 


oscy’*t OOSE*C oséc°e o000c*» OScl’sS 00S0°9 0S26°9 
0052°6 00F2°6 00S/ “6 oosl*F 00%/°6 0072 °6 00c2°6 
£999°6 SETb°S ¥BB6°8 6€se°3 O065°C vcoGs’9 90ZE°S 


OOMT*S6>T OOARPM* HHT OOTB°GET OODT°SET OOOS°OET OOvR SET OOBT*TZI 


OOBE* 86 009% °06 oor rs 00c6°E?L 0009°SY oorc 2s 0096 * ES o0v9 OL OOCE SE 
ooog"te 000c *&I 0008 °91 000¥ "vl 0000°<T 0009°6 000c°2 0008°> o00r-c 
QNOLNOD ANI T OFON 
dI1l 3SON WOAZ NI B2SS°49 1SZv-étt 0009°S9 £20¢° 91 
EYYNI SCB9°PELEE £060 °6CB> OLOS* LZYBL? BE60 *066€E 
c¥YNI BSbb 0022 v8sB°oczsl 39°56 °960¢6 TE99°EL01 
dil 350mN WO KI 6190 °OL OTGr° EST 0009°S9 BIOe’vT 
c¥¥NI T60T°Cobe coso°’Zér B66C°CC9! BBSO*cze 
7288°C L66E °C £990 °% BOeC*T 
NI 00OB*EST 0009°9> o00c °EB 0000°e 
JINOD aSGNITL9 3ayo0 
weeee WIL -2-- sem ANNA LAY mmm Hem ADORNS LNAD ---- --=—- ASON ---- 
NOILINI3ZRO AGOB JITISWWASI x¥ 
L1¥U7 TINI NOTLYANSTANOD ANOH TWidWIS 
c 39vd SNOILIV8NYF ANON JTISSIWN Sua SOUHIIW JIlwwNAOOaIY 
l 2S) ¥ 69/£ ATM ¥ WUILYO WZIISSIW OFlwwoLAy 3vSN 3HL 


SINIOd SNONNIINODEIO 3AdOTS SILYVDIGNI ¥ - Z10N 


0006°¢ ocrB’s 
00S.“ 00S2°6 
19G3°e OZeT°C 


¥000%°0 
¥00SZ2°6 
0062°G 
0000°0 1tOvas AOO8 


¥OOOC “EST 


oors°9TT OO9B"TTT yooor Zot 


0000*reS 
0000°0 SNOTIVLIS T¥NIONII DONO 


O1OaINZ9 AWNTON 


awnton 

VaAy AIIIIA 

OTOSINGIN vaay 

VAINY WHOINY Id 

Ollvgs SOININT 3 
HLONII 


ad VHS 


MISDATCOM Input/Output (Cont) 


Figure 24 


62 


(Cont) 


O35 


cé&ce°o 1Z00°0=- 169Z°O 0600°O6 30E0°O 000° 
TZTE°O 1T100°O- ctec’o 6400°O TECO°O 000°e 
CZOr°O $000°O0= Sv6E°O <c200°O ASTO°O 000°C 
Z£€6c° 1000°0- T¥tc°o 9000°0 9800°0 000° 1 
St6c° 1T000°0O- LOBT°O ¥000°0 c400°0 006°0 
1T9BC°O 0000°0 8S60°0 1000°0O BE00°O 000°0 
303 SNOIStA-wdD IVI LINILOd-WD SNOISIA-NI IVT INILOd -ND wHd T¥ 
Sr¥089°O = BOLJ¥32 ALI TYNOTISUdOdd DVE0 AOI SSOZD 
SNOTIVINITeD JDAOI TWtxv AML NI ORANIINI ION ST INSWTADNT OVAO JSYA AHL - JION 
$.000°0- 00000°0 9BIb0°O 08450°0 000° 
¥ 1000 °O0- 00000°0 90170°O 082S50°0 000°E 
90000° 0- 00000°9 9BIrO°O 08290°0 000°C 
<~0000°0- 00000°0 §B170°0 03250°0 000°T 
10000 °O0- 00000°0 iS 9b6TVO°O O#ZL50°O 00B8°0 
00000°0 00000°0 g93TRN'O 0U460°0 000°0 
VHA T¥-9 3ISva-v) JANVA/IENSSITEd-VA NOTIII 8A-¥9 wHd T¥ 
000°0 0Sl°s9 00S °6I 00S AT 9¥9° BGS 00°0 96°O0- 9N+IOET*& 02°90 
NI NI NI wT Z¥¥NT 130 430 13i/l &8 930 T¥¥NI /61 338/13 13 
TWILLIEIN "ONOT “Iv “OND VIAY J 19NV A TONY s3awnn s3guUnn 
S3INID “338 INIJWOW HIONIT ° 358 "338 1108 dF 1S307S$ SOIONAIS ISNIVSIdw31l IsENSSIId 41130130 J0NII Liv HIvwW 
woresennn---— SNOTSNIWIN TINAATIIE ---------- = wwe meen ween cee meee wee weee=s SNOT LIAINDD INDI 14 oon e nnn mw nw n ene eee e censor ee= 


INdinO WwIISvd JNO lV i008 
L¥7 TINA NOTIVANOTINOD AUOK AWdwIS ' 
t IDvd SNOI LYNN) JINOD JITSSTW ANZ SOOHIIK JluvnLoogssay 
t 3sv5 Y AR/L ASA ¥ WHNDIVO BVIISSIW AALYWOLNY WSN IH] 


MISDATCOM Input/Output 


Figure 25 


z coe’ él ATE*O ZOT°O c£0°O 00° 
066°0¢ ee” TOTO £20°O 00°t 
Ov6 "Te &E1°O 00T°O vIO°O 00°e 
90C°ET €90°0 OOT°O 900°0 oo°l 
99¢°Ce 1$0°0 OOT°O £00°0 oR’ Oo 
vec’ ve 000°90 O0T°O 000°0 00°0 
“d°I-x O2/15 a2 13 VWHdT¥ 
00+3000°0 TO*+FTABO° T= LO-TvHI°S- 10+7620°T TO-3PLO°" 000°0 660°0 G00°0 
00+3000°0 TO*3TEO0°T- TO-F30U°S- [O43 SHO°T 10-290£°9 000°0 £60°0 £00°O 
00+7000°0 TO*4¥60°l= l0-IBAG°H- 10+37T60°T T0-76039°¢ 000°0 960°0 v00°O 
00+3000°0 T0*3960°1T- TO-38C2°b- 10°3S60°1 T0-2060°SC 000°0 360°0 ¥00°O 
00+3000°0 10*3960°T= 10-3699" ¥- 1027960°T TO-TINL°Y 000°0 360°0 ¥00°O 
00+3000°0 TO*3TL60°1l- TO-399S"b= 10°3L60°T 10-3920" 000°0 960°0 b00°0 
$11) ANTI RAI vw vn) 119 NT) AQ 
TWNOTIIISIC WaILVYT TYNTONIT ONO eee TWNOTLIZEIC TwealYT -<-- 
coe ooooooorococco= (NYUIOYUN 83d) SINTLIYVAIRI omer em -eneerenrnrnnn- 
000°0 OS2°S9 00s él 00S °6I B¥9°86c 00°90 0f°0- SO+90ET°Y 
NI NI NI WI Z¥YNI 5an 930 Last 4 OAC 
WITLsan “SNOT “Iv “ONNT viv a TONY 2 TONY SIGWAN 
S3JINID “JI INIWOW HIONZ1 ° 378 “37a 11048 ALISAGIS .SMTIUNATE FJANIVAIIWAL 


SMOISNIWIO JININIIIE 


SNOILIGNOD JHOI13 


660°0 
0o1°’od 
oot°’o 
00T°o 
00T°O 
OOT°O 


v3 


S¥YNT/AT 


VANSSaad 


L920 
elZf°O 
TBO 
Tet°’o 
ESTO 
000°0 


wo 
IYNTONITSONOT 


e 


338/14 


ALLIOTVAA 


6£0°0 
£20°O 
£10°O0 
RO0'O 
900°0 
000°0 


NI 


i3 


AMNLITL Ty 


ae eae ewe e@eeeceeae es @e ee ee eee era @eoeeee" 2 


00° 
oo°e 
00°C 
00° T 
08°O 
00°0 


WHd1¥ 


0Z2°0 


sIGWNiw 
HIVW 


SILLISTSALIVAVHD JIWYVNANOZAY JIIVLS ANOTY ACaE 
1¥1 T1NA NOTIVENOTINOI AGHE ATdwIs 
SNNTIVANOTINOD ATISSUIW AQI SMOHLIW JINVNATOSIY 
¥ 68/2 ASS ¥ WODLYG STICSIW GAJVWOIMY JYSN JHZ 


¥ J9vd 
I asv9 


(Cont) 


MISDATCOM Input/Output 


Figure 26 


64 


10-3699ZT°I- £0+¢36¥S"*1 OOFIILZTCE o°y 
10-326650°1- CO<36fC°! 00+3Z98°C o°e 
ZO-III1val*sb- FOFTOILT: [ OO*7CnS °C os 
TO-3CCEIS°B- CO*3CCT 1 00+2660 °C o°t 
20-379266°L- £O+ILSO°T OO+A4HST°? -8°O 
70-3091 L4°L- 70437696°6 OO+3bVO'S 0°°0 
UYNI+OND NYND OND vid TY 
err re re nnn enn eens (NYT OYUN 83d) SIAIIVAISSG JIWWNAQ on---------- 2 --- =e 
000°0 0SzZ°S9 00S’ 6! 00S "61 8r9° REC 00°0 os*0- 9O+IT0CT°> . 02°0 
NI NI wl NI T¥¥NI "30 %30 last a 930 Z¥YNI/61 938/13 li 
WIITLdIA °9n07 "Iv! “9n01 vany AT9NY 119N¥ a3GHNnw SIGWNN 
S3IMN3D °338 INIWOW HIOWIT ° 338 "138 1108 4) 1S30IS SOIONAZS BENIVSTIWII BENSSASd ALIIDOIWA JaNITII¥ HOU 
moo eeenn---- SNOISNIWIM AININZIZE ------------- wn nn nn en ene mm n ewe n ene SNOTIININODD IH9IVS 2--- ener nen nnn nn nnn nna 


SJAIIVAI SAG JTWYNAG JNOTY AGO 

L¥VT TING NOTLIVANOTANOD ANNR AWdWIS 
$ 3Uvd SMOTIVANGT ANOS JTISSIW AU SGOHIIW JIWYNIGOAIY 
I IS¥9 ¥ 68/6 AT4 ¥ WODLYN JTISGUIW AFTLYWOINY AYSN FHI 


(Cont) 


MISDATCOM Input/Output 


Figure 27 


65 


LIST OF REFERENCES 


1. Mendenhall, Mee, Lesieutre, Bie Carugo, S Gar 
Dillinius, M.F.E., Kuhn, G.D., Aerodynamic Design of Pegasus, 
CP-493, paper presented at the AGARD Conference Proceedings 
#493, pages 7-1,2. 


2. Air Force Flight Dynamics Laboratory, AFWAL-TR-86-3091, 
Missile Datcom, Volume I- Final Report, Vukelich, S.R., Stoy, 
S.L., Bruns, K.A., Castillo, J.A., Moore, M.E., December 1988. 


3. Air Force Flight Dynamics Laboratory, Technical Report 
AFFDL-TR~72-83, Vol. I, 1972, Prediction of the Six—Degree-ome 
Freedom Store Separation Trajectories at Speeds up to the 
Critical Speed, Vol. I--Theoretical methods and Comparisons 
with Experiment, by Dillenius, M.F.E., Goodwin, F.K., and 
Nielsen, J.N., 1972. 


4. Joint Ordnance Commander Subgroup For Aircraft/Stores 
Compatibility, AOP-12, Volume 3, Aircraft Stores Interface 
Manual, pp. 2-(15-18) ,@iereal 19¢ee 


5. CricketGrapn” software, copyright 1990, Computer 
Associates. 


6. Nelson, R.C., Flight Stability and Automatic Control, page 
113, McGraw-Hill Book Co., 1989. 


7. Dickenson, B., Aircraft Stabiriity and Control for Pilome 
and Engineers, p. 14, Pitman Publishing Corporation, 1968. 


8. PC-MATLAB Version 3.2-PC software, copyright 1987, The 
MathWorks, Inc. 


Oe Oe software, copyright 786, Systems Control 
Technotogy,. Inc: 


10. Spahr, H. R., Theoretical Store Separation Analysis of a 
Prototype Store, Journal of Aircraft, Volume 12, Number 10, 
page 807, October 1975. 


11. van den Broek, G. J., The Use of a Panel Method in the 
Prediction of External Store Separation, Journal of Aircraft, 
Volume 21, Number 5, page 309, May 1984. 


12. Miranda, L. R., Extended Applications of the Vortex 


Lattice Method, Vortex Lattice Utilization, NASA SP-405, page 
27, l9s8ae 


66 


Mmomebertim, J. J., smith, M. L., Aerodynamics for Engineers, 
Prentice-Hall, Incorporated, 1989. 


67 


Oe 


INITIAL DISTRIBUTION LIST 
No. Copies 


Defense Technical Information Center Z 
Cameron Station 
Alexandria, VA 22304-6145 


Library,- Codems2Z Z 
Naval Postgraduate School 
Monterey, CA 93943-5002 


Chairman il 
Department of Aeronautics and Astronautics 

Naval Postgraduate School 

Monterey, CA 93943-5000 


Commander Ai 
Pacific Missile Test Center 
Point Mugu, CA 93042 


Naval Air Systems Command 1 
AIR-530E 
Washington, D.C 2036 


Flight Dynamics Laboratory 1 
Air Force Wright Aeronautical Laboratories 

Air Force Systems Command 

ATTN: Mr. William Blake 

Wright Patterson Air Force Base, OH 45433-6553 


Naval Weapons Center i 
Code 3592 

ATTN: Mr. L. Gleason 

China Lake, CA 93555 


Pacific Missile Test Center 4 
Code 9053 
Point Mugu, CA 93042 


Prof. 0. Biblarz, Code AA/Bi 1 
Department of Aeronautics and Astronautics 

Naval Postgraduate School 

Monterey, CA 93942-5000 


Prof. L. Schmidt, Code AA/Si i 
Department of Aeronautics and Astronautics 

Naval Postgraduate School 

Monterey, CA 93942-5000 


68 








rvatiier Mee O 
y rel , 
dat T ‘ i= 





Thesis 
H201235 Hansen 
Cal peo Sao Sep. otto 


methodology analysis. 





















































— "= s ei 7 5 L ty | 
offi! 4 bt ; r ; j & - a ae i 
L) t eo 4 th t rt 
aan ' aia * ef M oR 
a ig ha ar) - 
fd aa O . Ht | aa ar j 
é r Et , ¢ ry , br a A ‘ 
t,t. n F es , é 
r) UL “a ict tt hg a e F 
ay é Orta i ee oe | ae 7 ra 
ig ~ : ; Oh 4 : f] O ® is ar 
) j e, A )#) / } i , st | . a) ry 
} £/DO ee 2/2 VU FCM aA CAR Dn mk Le : 
, aie alt ry 44 whoy rw 4, Fa L. 
Tr “ity i a a: PP reg ' " 
e: y ahs P 
ue es PASO he As , mA 7a am ad 4 ais 4 ' : 
Lo el x at ie Ce tags os ‘ r 
2 PE eo é& 2¢ 8p ries ALY) L ‘oe. of a ’ 
as Te rs VOLS XS ONY i tn Pat Bie Sal a Sheds |! Tm 
Yop, ge ay are wy fer Le a a he ‘at oy ee Pe a 
OAR ey RAF d We tre uh A Py eM Say” Oe ee 6 ey fee ot AS aie e 
re fe, pe a AAD hy HO Ae e “(Gee Ar | oF ve 4 An P ee i. et vs AL ep , - 
pete re, pee ae es eS Pa pra s Bt ee Pre ACA ae te PaO ws me ve I Ms At Ie me ie Le a Ly 
TET wre A Pyar + : td ks SM Ai te Aa SY AA OO Hanh athe nelet t,o ae eu ee, y SF a ie 
Payot yyy earn PP eh: iy Rt a uf Orbe pe kiges 8545 44 fae ge og at "Wh Sige ty 1) "eee + P 
2 pe ye) ‘i , Lae wy FP as,° & \ ay ad ar Mar eee A oy oof 
es EO Pik Pye EAs BUG! A on ta Fi Cra ae ) A S in ee he ig me r hs ew a 0G 
Porter er ‘ be} talked on id dl hay CL are, Ply 7 eels hee? el Py A my ] Aion q ad ; . oe s i > rote £ re a ad 3 F . 
Sie artae Se oug NP ete cant, PO ae ay TLE ey ee wer 1 adits rw) he's , 7, f 
Very fe roe aie : rfl 4 a ee we ek Aas ; s . ‘4 s P ‘ F, PD a) rd 
iy) Auf ty, PE A ae rer en sy i a re ree. ee , 
en cee : c Pe AR RO CO a SU ar ad fence anf Pat ae Py 
enti rv, 5 - Po RO SY oe r He yi f Re CMMI. nL ot AUR Se Be | rary) Pi Fi : 
.y Pr ee ey VO EA Se Way fhanennt ‘ge r Che oe ° A , r 
oN Sy Warn A PKL ke en) a oY) is FALE OO Oty ery eS TM hs Peli sen he ged hd eet [you j 
ey) EATS MESES Sy a i hd , Py Sk ae ie AF 7 yr DA PAD ee AO a ee), SO fee eae a fe Fs Wl, FOP Be te ee a \ ~ , 
ape CW es A hed ME SS YS ee Ft SPE W 0 6608 Oe ogre wig oy of as YP be wy aU Plame Fim Sak de oe BL bed F P 
eR CT te} é ti ee Te at Poe PO r ¢, ay j ’ a Pr Ren 4 Z | | P Ce Par, S P 
r ees ee Bf Ps eve tees, yy oC Wa ay PTTL YO a aU PEP fe me 4 
ar Cee Pity Om ug een ONIN AU RD ay Pane ee) Baie Ml oe tel. A ee - Ry 
gh Tere pty ihe OR MDa AAS ERE OC RO ema ey Oe oe itp hte iets iO) aL he y 
aCe Y SAN eee s Se Ceased TPL ae) aaa PP Fk eS i Po) ‘ Cee gee en ae f P ae } 
ce be a a) : t Oe wae . 4 a af A ae rie ae P ] 5 
eas re oP: Cre, Oty ‘ ey : Pra ha ae F & x ra ON Le ae ee ‘we i 
Ben osath RE ote eaegh § re Be es ) ay 5 
is r Ae | $3 cs Py t Fa 
VF e a re Ls wes 4 oh a a Pe Uv : a 
. aaa ar a ace a rE, rf | é fee a e Eh 
Ra rt J a MIT Be aaa oo ; 
ra qi ate tote ree 5 Pay abe ar U ar a ( f 
; PoE mee j f Ais tt ‘ug. aoe 
J ‘ hea S & ol aT etme , 1 a | ee er es "4 5 ' 4 ys oa 
r) ‘eat 4 ome 5 A Fado +t , ; Pal ar a v's Le sore eee wae ve A oy , 
; TS ny, tee YEA A FL ORE irae PLT RON a) a Ben Pay. 0 a eh een a OT - ¢ 
eet ee i one xg pt FU AY POR AE ra 7 ey , rap PEAS: va) t PEPE RIA IG OST < 1 : oF sett? oe ie Ss ara et ea 
ee eA Ppte " : , me Pugh | ! Mr Se oe Oa pe BY ry rr i Pt oA ta at ace hy a eC A Pr Se ok ieee ie ri : 
ee ph ear one é * oe a p Be - a) "eat i wis 1 4 By * oe , Ce Pe ae ae i CSR ag Ly , at esate re ‘ 4-4 J FF ae he ‘ 2 , of rs ’ 
aes oy pipe A FD Ee 5 Mn de aes Lea te 5 ¥ 6 “gag ae >t + 5 re or, Lf © i hae " Pa Pa Aha 9 rd fel ‘af Pat LEP gh ah , ihe « ere Meee: i ie a u ' : ‘ 
re Fes wha i £ eB lt tie ste tits eg ity Acne é: , igi tetg 6 an? yf ve | a t 
RR Bee RA CRO a ee EMU AHEe botetha gies ANS g har eae © eS Pre é Fi 
o~ TP oe 1 ea eee f s UP ta Chea We athe rarer ey - 7 ; ; Hs sgen 
ae eee cr OL OP 2am! PUPS LaLa el ICR FA FSA ae J F ee Fi t 
Shee iron ta r} ion ots Wy, e's ; Cy ere AAW AAD F a , rn hr zi 
(i tebe ht A Ler t at BB 3 f : + = F ; "6 
reir wee ete om i bs ‘ Pah, ay SP UA Ee. DR ie aa § al ir a ? ; % ee 4 3 
/ ae A me ha Kf et Sat iis Irby be ‘ 0 5 AU 
oe ae, Bae 4 ee atte Bs fe ea ate ne A ee os PPS aT LSA SMES : i SUT ERP aT ted AOS DL at eR Bees 4, ar)! ae wey p 
Fee ’ F : B x BSS 3 fi Sue y : a ‘ 1 ; A 
rr y ¢ 
* pees 










rT) 3 Y R j "5 ; A aa 
* we rf L A ry 4 *i 
Sore ae ive nets Bana DoS Oak ie 4 ne 
é EE PSF RE SS ral eat Dat re LE ree LW ep? 
bet Sano tpn te ee Mehr n 
=. i 7 co 
a yah a as 
Sie ae = Cre ty oF i alates 


ae PET aie 


nee irre te iy 4 


Oe 









ea nT 
re a t ¥ 
3 lie a3 i €: 
< of rs Ly e fy tes heh ee 
a A rd ric a fs 
oan wld a IR Fy | 
Sh - FA *, 0 fs% Ts gs 
Ae a de eT 2) ru fe e . ry 
LY Fs { ¥ 5 
Ps r O) 
sf $4 








ee ae 
2 ee 


re er 
or 


Saar EN “2 
be ve at at Ue 


yay % a see eg! ° 
She — CaM a5 is , 






bs thd oda 

fy RO 

. AY ra 

AR Seb 
Lie Anthea hatte ett Petite) el Py ot Uf S9 i 

Petar) Sa Lott are eh vote 4 Shy 0 

ery i) 



























































































3 
fy 7 ‘ i ry 4s . On ir) r , 
“wed, " aT) a) 4 , " a es A Car 
i eS Oe 2 o¢ * PT} ers ny 4 & 
ar ahh: BRT wre ade i ye LEC we ae ei : hy 0 aos ‘ 
Meade eS a Le ee tf ¢ Pt 4‘, Ate é 4 ‘ 
A Sst ths €, A + a aod SS h ce un | " 
FEY $0 CaS ed ah eo We Seat 2A OG ads et ar TAO : te ea < 
ct 3 A Cm Pd geen 4 art ar} 
4 y ) i Se a es 75 ‘ ° a 
ates % mig eusnen a rh ‘ of or ve Le oe <4, , ry - 

7 o Vee Ge sis re Ph es é A 
apa eas 3 : ; LY 4 Sd a P é 5 5 ‘, 8 Ss in 
an a H, ng) mae 6 We i Ue ate ae: JC 

Ace TER: opie i Se of 
b hel eta ea had - A Fp Lh ee ‘ f - 7 
eres Praha ray mo y . ‘oe ae a , - P a4 ’ 
Med glans de Tans Le i } cn ara A re | rere 5 eS 7 we e t . 4 
hee 8 Pe ee a any y % © A a a Nae $2 ® viel os bo 7 is i A ° 
“Vit Rees ae Sy oF Ay we a: rf ata oF ; oC 4 o - ae Aes i ne a 
a Gey ph < LY epee Led retin ae Rrachih) cis ‘ eae ee R Ls ATA ore yey ¥ % _ ke Pili Fi Ph Pen he OM eatiy a Th OG parma 
penetra ot! Gah Bins: ct Say bred rye Sn ee ww yas in A au 5 re i wy at aL), P &e eee er a rt : r » 
Ph haa! wih wa ; ; He : we yy oh yh ads La tan i a Pate . ny as 5, ‘e. ree ‘ He e i” ts Fs > : fh: a ee eae , *e 
O™M 4 Me vee Se Te 4 , ae t Mi : teh, A ° a Eh r A ; ‘ 
oak ny \ CS ae are SSIS Reh . eh Ker a a ae a i i n't “4" ya WA Paraaa . hs on f ; ae oP F , Pee R 3 ° 
Ele! + ited Se: net hd Dee Shh ot \ Y 4 7 re i oe Ty ° 4 
aie ear! whe eth ba ih eran 5 5 [Ss A ‘ es Hd he iY 45) wei: Op Paar vat OR NA Np & * iy et ye i Aes a $f P Fie ee A AY ri aa O Ei : Fs ny é 
f a Zs ep bs j ri - . 7 ‘ as, - n ; 
a ua eS >! ay nN ie a i. 2 4 sor} ban , a Tate te Po va ore Rs ag yt at sn 2 aa ey ra be < Fe i) Si “*, oats cs Tks eee 5 ee oi yo , 
by res 4 = y , ark Ae a Mats St ree q ey E a y F as : 
aS asites uno sey pei ven eet ad . , A  ehae NS 4 br Cae Daren ge Hae 7 ra . vs * P , a A as ai rh ee ree an) 
ho Letty be wy Si hehe nena bars abe A hg \ - DOE Pre BN ha bgt ra ice o Cleek ~ ie Cae 8 S ee Be Pe 
i ; cae Pg Dey Daten et ‘ a 4 Js) Ld +e 6.U Agra, & a Phra’ wn MT af) Re, aa F ae P 
; eee Le tk Fetes BAR ! pe 3 Shee ee Nn Feet Re A aatel? ible as AS 2 ‘ SN Ae SF Fi Sa ott ee Ls wh fi ; 
I oa te A o ¥ nt D 4 on) © ay ry a ia z m o 
DREN yaaa NH ae eo eR SAP Clee ari 
Cae Re et ¢ SER Por par aes 4 ae WO ay , ° d . A $e 
fa? Panay Vedi taad | Spy ie as a 7 wy, a PA ee tet A) ded Fe Sa ey Te ey Ft al er a “Tae f f é 
ts eet pb th Meet vay 4? aR Tet bo ipa ea AS ek Rare Ons OMe nly og’ dE OO Oey en me i) a 
| S 7% a | \ ® s 
potetae RK x eee 4 vane ft es Be eS RCA ris WIA alas tal Soy On i ere gE re j Re el oar F * 4 ca 
royal: yy ere e ae, ooh secre pay atte Aen ein db. bade, ta 0 en TY ries Ra DE, Ae oY A a eee r i a e@ x a fi or 
tap ams tt FP Co A :¢,8 ig AH hs wah tte Lat y mn J Fb ye aA i A 7) ’ | a “ar ee ts : ay ST is ra * P 
, Ee =) b AS th Re wee ops auth Pp ae aa ae Pa PO ‘ Pa: ' 
wees ads. ated ey ry Pato Pale , os ze ee eh Co Be Dg 4 " tn a PhS 5 
Lay Sy, Yah Xgl hott et Lod dad Jat tek te Bot pay Vian 5 teh is a ad a - “ven Rr ie 5 f +4 PL eh ‘ ae , 5 
Racks AA es atte Ae ’ nes i CT Pa tha ta CS fy aria h £ ar edits tt e BEG Le a) ° 
bendy Shitethen 41k Y ap eae alata ed i 9 bf ; Pde) eS bd ‘ie! 610.¢ b.wt & so f +s, 1 ey st: - Po 
“ag > as Pathe Fo La + € Je A) mB* @me @ b i : 
p RETR : Oe ay tial A R fo & eae See at ; & 
Papillon, “an . ti ny rena ” sue Are ‘ o bu Be e ; E peeve aS E ot a Pee - bf 
Neu eh yee t C Ek oh ie a ed ee wear CUE Bs oe ; ‘J =. CY , i F 
By 7 i A | 
eis Fath POC vi i Lm} fig Py a ay ; ar - ari ° 8 
=. i + eh hy a ae | ee err = “° 7 C2 - t Ee y P s 
tat a oh ee ELS ee a} 3 4 a) Wousale bd ba bd b . Ly tr} ry 
Fae CD te Se Me ek Oe ee 9 a ” m i", in} i Li as ry 
eee x reas ef Pg ah i ee ih reeled | s0ye et ees 8 = hi aM a a's © a ee | ete iY 
~e fe ae . 7 Fl ° “aes a, a Py ny 
Se Sw ie ear ha Pt ae ar ae kT hy 4 =e om a. 7 ‘se8 ns 
Ra ee ae - ae ATO Ci 28 ar ata “ 5 " a? 
©. via Se hs Carey oy Ate afer 7 1 . * Py . 5 A Cee ae Pe C 
Lda he! 5 Chea i rr T fi A - oth, ee 3s a! : ‘ ‘ 
bit ek ath et ee P os Por fe s r r : ¥ 
he de Ct ee coe thy Sate E é , . a ss HY 
Salts ot Ma ti 7 aw 5 e D A v ‘ ‘ Cr rs mn +4 
ot , Px oe On ee aoe ge eC ar F 
othe Aa &% re ear a 5 P - rn A 
at Yer . a acer ote, ~ 
Rediv* saapesal AS ee ete 7 is . i de bt he M bd 
Ss Sue a : PE SRA ys ty roared an se ioe . 6 Om ie Pt s se a i th Fok rs 
F : ra 4 Ti a ere wtat a? i A> a Pw * Ti ae ea n P od . 
Sate oy TL er R PCO Ser An a lal SM re a Se ek: AA eT : o. 
A Care - F ' a 7 SCS ‘ 
gt LA AL x LY L- AA ty < hy yt vs . bd ) ca b . oa : Ls 5 ny 
Fad iM . a yr 2 < “fh 7 s 7 5 A r € fi ° P PF ms rs 
Vac Get oS H z 4 ig Am) Hi 
aS RPT ee es 8 fae ge ae TES tae es peta 4S ¢ i - 
rn ove eR te z . : Pa . ae | ° ¢ r Ps id 
a *ete A © sad? , why > eed a ¥ e Oa v z 5 Li ¢ Py 
ies "y > Oye ho Lace ys A oe a rt, Se | i Pp a e 
high ied aA, a v ie , ht, Hohe Rat Soa Le D .” r . 
re fe TH ret | tie * % S 1° st n hor rn 5 
ree Tg bar be ties cs a ry - P : 
et WOR Oo “oe AR ee Ae ew ie a 8 ; “ Py - F < F 
¥en. > 8° ee ee Te he a rar - 


